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Abdur Rehman 



Abstract 



We study the spectroscopy of the states which defy conventional cc charmonium and bb 
bottomonium interpretation respectively, and are termed as exotic states. In August 2003 
a state A(3872), was discovered [K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 
91, 262001 (2003) |arXiv:hep-ex/0309032 |] and in December 2007 a state F(10890), was 
reported [K. F. Chen et. al. [Belle Collaboration], Phys. Rev. Lett. 100, 112001 (2008) 
[arXiv : 071 . 25 77] . One possible interpretation of such exotic state is that they are tetraquark 
(diquark-antiquark) states. Applying knowledge of non-relativistic constituent quark model, 
we calculate the spectrum of hidden charm and bottom exotic mesons within diquark- 
antidiquark model. We investigate that A(3872) is 1+"'' state of the kind X\^cq] = ([eg] 5=1 [eg] 5=0)5- 
and F(10890) is 1 state of the kind y[bg] = ([feg]s=o[^^]s=o)p-wavc and calculate the de- 
cay modes of these exotic states, further supporting X(3872) and F(10890) as tetraquark 
(diquark-antidiquark) states and resolve the puzzling features of the data. We study the 
radiative decays of these states, using the idea of Vector Meson Dominance (VMD), which 
we hope will increase an insight about these tetraquark states. 



Chapter 1 
Introduction 



In Heavy Quark physics one of the interesting problems is to determine the properties of 
newly discovered particles. It becomes very appealing when these particles do not fit into 
the existing models. Quark model provides a convenient framework in the classification of 
hadrons. Most experimentally observed hadronic states fit in it nicely. The states which 
are beyond the quark model are termed as exotic. To interpret these particles many new 
hypothesis are created. Non-quark model mesons include: 

1. exotic mesons, which defy conventional qq interpretation; 

2. glueballs or gluonium, which have no valence quarks at all; 

3. tetraquarks, there are two configurations in tetraquark picture, molecular models 
and Diquark-antiquark model; 

4. hybrid mesons, which contain a valence quark- ant iquark pair and one or more gluons. 
The exotic states in the charmonium spectrum have been found experimently, few of 

them are labelled as X, Y and Z states. In 2003 a particle temporarily called X(3872), was 
discovered by the Belle experiment in Japan [1]. These exotic states refuse to obey conven- 
tional cc charmonium interpretation [21 [3]. So, many theories came: molecular models [HE], 
more general 4-quark interpretations including a diquark-antiquark model [6]- [10], hybrid 
models [TT] etc. To explain the nature of the X (3872) , it was suggested as a tetraquark 
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candidate. The name X is a temporary name, indicating that there are still some questions 
about its properties which need to be tested. The number in the parenthesis is the mass 
of the particle in MeV. There are two configurations in tetraquark picture. In the first 
configuration, binding each quark to an antiquark and allowing interaction between 

the two color neutral pairs [q'a^"] [q'/?^''] • This is what we call the molecular model. In par- 
ticular the X (3872) happens to have a mass very close to the DD* threshold. The binding 
energy left for the X (3872) is consistent with E ~ 0.25± 0.40 MeV, thus making this state 
very large in size: order of ten times bigger than the typical range of strong interactions. 
These questions apply to other near-to-threshold hypothetical molecules and have induced 
thinking to find alternative explanations for the X (3872) and its relatives. In the second 
configuration, binding the two quarks in a colored configuration called diquark [qq]^, with 
antidiquark [qq]°'. This configuration is what is called diquark-antiquark model where 
diquark is a fundamental object. The X(3872) is a [eg] 5=1 [eg] 5=0 tetraquark. The work of 
this dissertation is to discuss the lowest lying exotic meson X (3872) and higher mass exotic 
state yb(10890) in the frame work of diquark-antidiquark model. 

In the diquark-antidiquark model the mass spectra are computed as in the non-relativistic 
constituent quark model. In the constituent quark model hadron masses are described by 
an effective Hamiltonian that takes as input the constituent quark masses and the spin-spin 
couplings between quarks. By extending this approach to diquark-antidiquark bound states 
it is possible to predict tetraquark mass spectra. The mass spectrum of tetraquarks [gg] [gg] 
with q = u, d and q = c, b neutral states can be described in terms of the constituent diquark 
masses, mg, spin-spin interactions inside the single diquark, spin-spin interaction between 
quark and antiquark belonging to two diquarks, spin-orbit, and purely orbital term [6]. 

In the second chapter, we briefiy discuss the Quark Model, hadron spectroscopy and the 
concept of tetraquark. As diquark is the fundamental object in the diquark-antidiquark 
model, a complete section is devoted to understand its characteristics, especially parity, 
color etc. 
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In chapter 3, we give a formulism of diquark-antiquark model. The exotic state X (3872) 
is the focus of study in this chapter. We calculate the spectrum of hidden charm states using 
this model, which automatically shows that the 1^^ state is the X (3872) . In the last section 
we discuss the concept of isospin symmetry breaking which helps to understand the finer 
structure of the X (3872) and finally we calculate the decay widths of X (3872). We calculate 
the radiative decay widths of X (3872) by exploiting the idea of Vector Meson Dominance 
(VMD). 

In December 2007, the Belle collaboration working at the KEKB e~^e~ collider in Tsukuba, 
Japan, reported the first observation of the processes e~^e~ Y^bg-] — )■ T(15', 25") tt+tt" near 
the peak of the T{5S) resonance at the center-of-mass energy of about 10.87 GeV In 
the conventional Quarkonium theory, there is no place for such a nearby additional bb reso- 
nance having the quantum numbers of T(5S'). An important issue is whether the puzzling 
events seen by Belle stem from the decays of the T{5S), or from another particle Yb having a 
mass close enough to the mass of the T(5S'). The puzzling features of these data are that, if 
interpreted in terms of the processes e~^e~ — > T{5S) — t- T(IS') tt+tt^, T(2S') tt+tt^, the rates 
are anomalously larger than the expectations from scaling the comparable T{AS) decays to 
those of T(5S'). 

In the last chapter, we modify the formulism of diquark-antiquark model to calculate the 
spectrum of hidden bottom states for L = 1. We are able to show that Yb is J^'~^ = 1 
state, with Y[bq] = {[bq]s=o\b(l]s=o)p-wa.ve, with the value My^ ^ (for q = u,d) equal to 10890 
MeV We identify this with the mass of the Yb from Belle [13], apart from the T{5S) and 
T{6S) resonances. We calculated the leptonic, hadronic and radiative decay widths of the 
Yb that may solve the puzzling features of the data. 
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Chapter 2 

THE QUARK MODEL AND 
BEYOND 

2.1 Quarks 
2.1.1 An Overview 

A quark is an elementary particle and a fundamental constituent of matter. Quarks com- 
bine to form composite particles called hadrons. Due to a phenomenon known as color 
confinement, quarks are never found in isolation; they can only be found within hadrons. 
For this reason, much of what is known about quarks has been drawn from observations 
of the hadrons themselves [H]. Quarks possess a property called color charge. There are 
three types of color charge, arbitrarily labeled blue, green, and red. Each of them is comple- 
mented by an anticolor — antiblue, antigreen, and antired. Every quark carries a color, while 
every antiquark carries an anticolor. The theory that describes strong interactions is called 
Quantum Chromodynamics (QCD). 

The quarks which determine the quantum numbers of hadrons are called valence quarks; 
apart from these, any hadron may contain an indefinite number of virtual (or sea) quarks, 

antiquarks, and gluons which do not influence its quantum numbers [15]. There are two 

4 



families of hadrons: baryons, with three valence quarks, and mesons, with a valence quark 
and an antiquark. The existence of "exotic" hadrons with more valence quarks, such as 
tetraquarks (qqqq) and pentaquarks (qqqqq), have been conjectured but not proven [16]. 
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Figure 2.1: Six of the particles in the Standard Model are quarks. Each of the first three 
columns forms a generation of matter. 

In modern particle physics, local gauge symmetries — a kind of symmetry group — determine 
interactions between particles. Color SU(3) (commonly abbreviated to SUc(3)) is the gauge 
symmetry that generated by three color charges which a quark carry and is the defining 
symmetry for QCD. The requirement that SUc(3) should be local, i.e its transformations be 
allowed to vary with space and time — determines the properties of the strong interaction, in 
particular the existence of eight gluons to act as its force carriers |17] . 

Properties 

The Table 2.1 summarizes the key properties of the six quarks. Flavor quantum numbers 
( isospin (Jx), charmness (C), strangeness (5*, not to be confused with spin), topness (T), 
and bottomness {B)) are assigned to certain quark flavors, and denote qualities of quark- 
based systems and hadrons. The baryon number {B) is +1/3 for all quarks, as baryons 
are made of three quarks. For antiquarks, the electric charge {Q) and all flavor quantum 
numbers(i?, J^, C, 5, T, and B) are of opposite sign. Quarks are strongly interacting fermions 
with half integer spin. Quarks have positive intrinsic parity, and are spin 1/2 particles. There 



are additive flavor quantum numbers for tliree generations of quarks. Antiquarks have the 
opposite flavor sign. The charge Q and these quantum numbers are related through the 
Gell-Mann-Nishijima relation: 



Q = h + 



B+S+C+B+T 



(2.1) 



where B = 1/3, is the baryon number for each quark |18j . 

Table 2.1: Properties of quarks 
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2.2 Quark Model 

The quark model is a classification scheme for hadrons in term of their valence quarks. 
2.2.1 Hadrons 

According to the quark model, the properties of hadrons are primarily determined by their 
valence quarks. Although quarks also carry color charge, hadrons must have zero color charge 
because of a phenomena called color confinement. That is, hadrons must be "colorless" or 
"white". There are two ways to accomplish this: three quarks of different colors, or a quark 
of one color and an antiquark carrying the corresponding anticolor. Hadrons based on the 
former are called baryons (half odd integer spin), and those based on the latter are called 



mesons (integer spin). This is possible because of the remarkable property that a singlet 
exists in3(8)3(g)3 = l©8©8©10as well as in 3 (g) 3 = 1 © 8. All hadrons are labelled by 
quantum numbers. One set comes from the Poincare symmetry- J'^*", where J, P and C 
stand for the total angular momentum, Parity, and Charge-symmetry respectively. Hadrons 
with same (lowest lying as well as higher mass states) are distinguished from each other 
by some internal quantum numbers. These are flavor quantum numbers such as the isospin, 
strangeness, charm, and so on. All quarks carry an additive, conserved quantum number 
called a baryon number, which is +1/3 for quarks and -1/3 for antiquarks. This means that 
baryons (groups of three quarks) have B = 1 while mesons have B = 0. Hadrons have excited 
states known as resonances. Each ground state hadron may have several excited states. 
Resonances decay extreme quickly (within about 10"^"^ seconds) via the strong nuclear force 

m- 



2.2.2 Baryons 




Figure 2.2: Combinations of three u,d or s quarks forming baryons with a spin 1/2 form the 
uds baryon octet. 



2.2.3 Exotic Baryons 

Exotic baryons are hypothetical composite particles which are bound states of 3 quarks and 

additional elementary particles. The additional particles may include quarks, antiquarks 

or gluons. One such exotic baryon is the pentaquark, which consists of four quarks and an 
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Figure 2.3: Combinations of three u,d or s quarks forming baryons with a spin 3/2 from the 
uds baryon decuplet. 

antiquark (13= 1), but their existence is not generally accepted. Theoretically, heptaquarks (5 
quarks, 2 antiquarks), nonaquarks (6 quarks, 3 antiquarks), etc. could also exist. Another 
exotic baryon which consists only of quarks is the H dibaryon, which consists of two up 
quarks, two down quarks and two strange quarks. Unlike the pentaquark, this particle 
might be long lived or even stable. There have been unconfirmed claims of detections of 
pentaquarks and dibaryons [19j . 



2.2.4 Mesons 

The main difference between mesons and baryons is that mesons are bosons (which obey 
Bose- Einstein statistics) while baryons are fermions (which obey Fermi-Dirac statistics). 
Since mesons are composed of quarks, they participate in both the weak and strong inter- 
actions. Mesons with net electric charge also participate in the electromagnetic interaction. 
They are classified according to their quark content, total angular momentum, parity, and 
various other properties such as C— parity and G— parity. They are also typically less mas- 
sive than baryons, meaning that they are more easily produced in experiments, and exhibit 
higher energy phenomena sooner than baryons would. For example, the charm quark was 
first seen in the J/Psi meson (J/^) in 1974, and the bottom quark in the upsilon meson (T) 



in 1977 [20]. 



Classification Of Mesons 

The mesons are classified in J^^ multiplets. The allowed quantum numbers for L smaller 
than 3 are given in table. Interestingly, for J smaller than 3, all allowed J^*^ except 2 
have been observed [16]. 



Table 2.2: Classification of mesons 
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Particle physicists are most interested in mesons with no orbital angular momentum 
(L = 0), therefore the two groups of mesons most studied are the S' = 1; L = and S* = 0; 
L = 0, which corresponds to J = 1 and J = 0, although they are not the only ones. It is also 
possible to obtain J = 1 particles from 5 = and L = 1. How to distinguish between the 
5 = 1, L = and S* = 0, L = 1 mesons is an active area of research in meson spectroscopy. 

For lighter up, down, and strange quarks the suitable mathematical group is S'f/(3). The 
quarks lie in the fundamental representation, 3 (called the triplet) of flavor S'f/(3). The 
antiquarks lie in the complex conjugate representation 3. The nine states (nonet) made 
out of a pair can be decomposed into the trivial representation, 1 (called the singlet), and 
the adjoint representation, 8 (called the octet). The notation for this decomposition is 
3 (g) 3 = 8 © 1. There are generalizations to larger number of flavors. Charm quark is icluded 
by extending 5f/(3) to SU{A) [I6]. 

Types Of Mesons 

The rules for classification are presented below, in Table 2.3 for simplicity [2T] . 



Table 2.3: Types of mesons 
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Figure 2.4: Combinations of one u, d or s quarks and one u, d oi s antiquark in = 
configuration form a nonet. 

Flavorless mesons are mesons made of quark and antiquark of the same flavor (all their 

flavor quantum numbers are zero. Flavorful mesons are mesons made of pair of quark and 
antiquarks of different flavors. 

2.3 Mesons And Symmetries 

Most of the symmetries in elementary-particle physics are continuous. A typical example is 

the symmetry generated by rotations around an axis, where the angle of rotation can assume 

any value between and 27r. In addition to continuous symmetries, there arc also discrete 

symmetries, for which the possible states assume discrete values classified with the help of 

a few integers. In elementary-particle physics there are three discrete symmetries of basic 

importance: parity, charge conjugation and time-reversal. 

Spin (quantum number S) is a vector quantity that represents the "intrinsic" angular 

momentum of a particle. Since quarks are fermi particles of spin 1/2 {S — 1/2) and mesons 

10 



Figure 2.5: Combinations of one u, d or s quarks and one u, d or s antiquark in J — 1 
configuration form a nonet. 

are made of one quark and one antiquark, they can be found in triplets and singlets spin 
states. The orbital angular momentum (quantum number L), that comes in increments of 
IH, represent the angular moment due to quarks orbiting around each other. The total 
angular momentum (quantum number J) of a particle is the combination of intrinsic angular 
momentum (spin) and orbital angular momentum. It can take any value from J — \L — S\ 
to J = |L + 5*1, in increments of 1. 

2.3.1 C- and P- parity, and Isospin 

An important property of any quark state is the behavior of its wave-funtion under certian 
transformations. One important transformation is the replacement of particle with anti- 
particle, called C-cojugation. An other one, called P-transformation, is the switching the 
signs of all coordinates. Many states, but not all, are eigenstates of these two transformations. 
This means that: 

The former is possible for the states which are flavor neutral, e.g. electrically neutral. These 

numbers, Xc and Xp, are called the C-parity and P- parity of the particular quark state 
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respectively. These transformations have one special property: 



= = $ (2.2) 

p2$ = ^2^$ = $ (^2.3) 

It follows that the particular quark state may have C or P— parity either +1 or —1. For 
example, the P— parity of tt— mesons is —1: 

P7r° = -7r°, P7r+ = -7r+, Pn' = 

Here tt stands for the vr— meson wavefuntion. The C— parity for 7r° is +1, Ctt^ = +7r'^; but 
7r+ and tt~ do not have definite C— parity, i.e. Ctt^ = +7r^, Cvr^ = +7r+. 

For the system of a quark and an antiquark (gg)/, [note that q and q have opposite 
intrinsic parities], one has 

P = (-1) (-1)^ , C = (-1)^+^ (2.4) 

For states composite of integer spin bosons these formulae are different. For example, for 
TT+TT^ system where pions have zero spin and same intrinsic parities. 

P = C=(-1)^ (2.5) 

An important property of C and P— parities is that they are conserved in the strong and 
electromagnetic interactions [22]. A generalization to C— parity is G— parity G = (— l)"^^*^^^ 
for mesons. The isospin of u and d quarks is equal to 1/2, with u quark having positive 
isospin projection I z = 1/2, and d quark with I z = —1/2. All other quarks have zero 
isospin and same holds for a state which is made up of these quarks. For example the isospin 
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of any charmonium state is zero. In particular, isospin becomes crucial when we discuss the 
possible assignment for the X (3872) , which is assumed to tetraquark state. Isospin / is 
another property of quark-antiquark system, which is conserved in strong interactions and 
follow the same algebraic rules as the regular spin S. Hadron with nearly same mass can be 
put into isospin multiplets: 



Table 2.4: Isospin Multiplets 




Isospin is conserved in strong interaction and as such in that limit, member of each multiplet 
will have the same mass. The small difference then arises due to electromagnetic interaction, 
which still conserves Iz, since Q — -\- isoscalar and charge conservation then implies Iz 
conservtion and/or m„ ^ ma- This is illustrated by the following decay. 

eV = 1, 4 = 0) ^ A7I/ = 0, /3 = 0) 

Isospin and its third component are not conserved in the weak interaction, as demonstrated 
in the decay: 

2.3.2 Isospin Symmetry 

Isospin symmetry, which is an exact symmetry as for as strong interaction is concerned, 

is broken by electromagnetic interaction and/or m„ 7^ m^. Thus QCD Lagrangian has 

isospin symmetry. If m„ = = the QCD Lagrangian has chiral symmetry. Since 

rriu ~ nid « ^qcd, the symmetry of the QCD lagrangian is broken when = 7^ 0. 
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2.3.3 Isospin, Charge and Flavor Quantum Numbers 

The pion particle had three "charged states", it was said to be of isospin 1 = 1. Its "charged 
states" 7r+, 7r°, and tt^, corresponded to the isospin projections 1^ = +1, 1^ = 0, and = — 1 
respectively. Another example is the "rho particle" . Isospin projections were related to the 
up and down quark content of particles by the relation 

^z = ^[{nu-nu) - {rid-n^)] (2.6) 

where the n's are the number of up and down quarks and antiquarks. 

It was noted that charge {Q) was related to the isospin projection (/ z), the baryon 
number (B) and flavor quantum numbers (5*, C, B, T) by the Gell-Mann-Nishijima formula. 
Strangeness flavor quantum number S (not to be confused with spin). Flavor quantum 
numbers of composites are related to the number of strange, charm, bottom, and top quarks 
and antiquark according to the relations: 

5" = -(ris-ns) 
C = +{nc-nc) 
B = -{rib-ni) 
T = +{nt- ni) 

This implies that the Gell-Mann-Nishijima formula is equivalent to the expression of charge 
in terms of quark content [23] . 

2 1 

Q = o [('^w ~ + i^c - ric) + {ut - rit)] - - [{rid - ^j) + (^s - rig) + [rii, - ni)] (2.7) 
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2.4 Exotic Meson 



From the quantum numbers in Table 2.1-2.3, there are several combinations which are miss- 
ing: 

0+-,0"",r+ and2+- 

These are not possible for simple qq systems and are known as "exotic" states. Beyond the 
simple quark model picture of mesons, there are different frameworks suggested to accomo- 
date these states with exotic quantum numbers. Non-quark model mesons include: 

1. glueballs or gluonium, which have no valence quarks at all; 

2. tetraquarks, which have two valence quark-antiquark pairs; and 

3. hybrid mesons, which contain a valence quark-antiquark pair and one or more gluons. 
All of these can be classfied as mesons, because they carry zero baryon number. Of these, 

glueballs must be flavor singlets; that is, have zero isospin, strangeness, charm, bottomness, 
and topness. Like all particle states, they are specified by the quantum numbers J^^ and 
by the mass. One also specifies the isospin / of the meson. It is an old idea that the light 
scalar mesons a(980) and f(980) may be 4-quark bound states. The idea was more or less 
accepted in the mid-seventies but then it losts momentum, due to contradictory results. If 
the lightest scalar mesons are diquark-antidiquark composites as shown below, it is natural 
to consider analogous states with one or more heavy constituents, to be discussed in the 
following section. 

2.5 Diquarks: An Introduction 

The notion of the diquark usually means the system of two rather tightly bounded quarks 
with a small size of 0.1 — 0.3 Fermi [25]. This section is devoted to diquarks and their role 
in understanding exotics in QCD. Diquarks are not new, they are almost as old as QCD. 
Gell-Mann mentions it prominently in his first paper on quarks in 1964 [26j . 
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usus 

• m ^(980).i^j(980) 

dsus 

^ k(800) 

udsd 



9 ^(600) 

iidud 

Figure 2.6: Identities and classification of possible tetraquark mesons. First horizontal line 
denotes 1 = states, second, 1 = 1/2 and third one, 1 = 1. The vertical axis is the mass. 

Baryons can be constructed from quarks by using the combinations (qqq), (qqqqq) etc, 
while mesons are made out of (qq), (qqqq), etc. The lowest baryon configuration (qqq) 
gives just the representations 1, 8, and 10 that have been observed, while the lowest meson 
configuration [qq) similarly gives just 1 and 8. 

The constituents of the tetraquarks, diquarks and antidiquarks, have well-defined prop- 
erties, characterized by their color and electromagnetic charges, spin and flavor quantum 
numbers. The tetraquark hadrons are singlets in color (pictured white), and hence they 
participate as physical states in scattering and decay processes. This is not too dissimilar a 
situation from the well-known mesons, which are color singlet (white) bound states of the 
confined colored quarks and antiquarks. 

We follow the suggestion by Jaffe and Wilczek (1977) of having diquark as building blocks 

|15] . Diquark correlations in hadrons suggest qualitative explanations for many of the puzzles 

of exotic hadron spectroscopy. Operators that will create a diquark of any (integer) spin and 

parity can be constructed from two quark fields and insertions of the covariant derivative. We 

are interested in potentially low energy configurations, so we omit the derivatives. There are 

eight distinct diquark multiplets (in color x fiavor x spin) , which are enumerated by R. L. Jaffe. 

Since each quark is a color triplet, the pair can form a color 3c, which is antisymmetric, or 

6c, which is symmetric. The same is true in SU(3)-fiavor. The constructions look more 

familiar if we represent one of the quarks by the charge conjugate field: qq — )■ qcq, where 
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= —iq'^cr'^'y!^. Then the classification of diquark bilinears is analogous to the classification 
of qq bilinears. There are only two favored configurations. The most attractive channel 
in QCD seems to be the color antitriplet, flavor antisymmetric (which is the 3/ for three 
light flavors), spin singlet with even parity: [qq]^'"^^'^'^ . This channel is favored by one gluon 
exchange and by instanton interactions [2^ . It will play the central role in the exotic drama. 

\{qq}3ciA)3f{A)0+{A)) good diquarks (2.8) 
\{qq}3c{A)6f{S)l^{S)) bad diquarks (2.9) 

Both of these configurations are important in spectroscopy. Now we can construct the 
operators for "good" scalar diquarks and "bad" vector diquarks [0]. Heavy-light diquarks 
can be the building blocks of a rich spectrum of states which can accommodate some of the 
newly observed charmonium-like resonances not fitting a pure cc assignment. 

Qia = eijkeabc{i(T2)q^\'"' = eijkeabciqi^q'"') (2.10) 
Qo'' = eabM'lQ''' + tl<in (2.11) 

Both represent positive parity, O"*" and l"*", states. We work out the diquark masses in quark 
model, in which all residual quark interactions are incorporated. Diquarks are, of course, 
colored states, and therefore not physical. The good (scalar) and bad (vector) diquarks 
configurations are our main interest. 



2.6 Te t r aquar k 

A tetraquark is a hypothetical meson composed of four valence quarks. In principle, a 

tetraquark state may be allowed in Quantum chromodynamics, the modern theory of strong 

interactions. Examining the color algebra of the system reveals there are two independent 

tetraquark singlet states: 3(g)3(8)3(8)3 = 27©10©10©8©8©8©8©l©l. But they can 
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be obtained in four different ways, depending on the intermediate color states: the singlet 
scheme (or molecule), the octet scheme, the triplet scheme and the sextet scheme. 




8 ® 8 = 27®10®10®8®8®1 



1 / ® \ 1 /= 1 (molecule) 




6^® 6\=27®8®1 



Figure 2.7: Singlets can be obtained via four different schemes (only two are independent): 
octet, singlet (or molecule), sextet and triplet. 

A diquark is either in symmetric color state 6c or color antisymmetric state 3c. The 
antidiquark is either in antisymmetric color state 6c or color symmetric state 3c. Now 



6c ® 6c 



3c ® 3c 



3c ® 6c 



35c © Ic 



© Ic 



10c © 8c 



6c® 3c 



10c © 8c 



Hence only 6c ® 6c and 3c © 3c give color singlet state. 

Now diquark are either in symmetric or antisymmetric in flavor: 

[qq] = -^{mj - QjQi) h j = u, d, s, c 

{<i<i} = -^{mj + mi) 

For antisymmetric color state 3c or 3c, Pauli principle requires, overall wave function of 
diquark or antidiquark to be symmetric in flavor, in space and spin: (s, denote the spin of 
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diquark or antidiquark) 





; P = 


1 




; P = 


{(!Q}l=o, s=i 


; P = 


1 


{qq}L=i, s=o 


; P = 


Ml=0, s=0 


; P = 


1 




; P = 


{m}L=0, s=l 


; P = 


1 


{m}L=i, s=o 


; P = 



We have a nonet of low lying scalar mesons 0"'', composite of tetraquark viz 

[gg]L=o, s=o Ml=o, s=o 

As an example, we consider the following two tetraquark quarkonium states, with L = 0: 

=0, s=0 / 

The two states have J^'" = 1^^ and 1+^. For q = c the J^'" = l"*""*" meson is identified 
with the state X(3872). This state was reported by the Belle in 2003. It was suggested as a 
tetraquark candidate. This state was discovered in the J/\l' tt+tt^ distribution. 

In 2004 the D(2632) state, seen in the SELEX experiment, was suggested as a possible 
tetraquark candidate with quark contents [c(i][(is]. In 2009 Fermilab announced that they 
have discovered a particle temporarily called y(4140), which may also be a tetraquark with 
quark contents [cs][cs] [23]. 
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2.7 QCD, Quarkonia and Hadron Spectroscopy 



2.7.1 The Important Physical Properties of QCD 
Gluons 

The gluons, being mediators of strong interaction between quarks, are vector particles and 
carry color; both of these properties are supported by hadron spectroscopy. 

Color Confinement 

Confinement which implies that potential energy between color charges increases linearly at 
large distances so that only color singlet states exist, a property not yet established but find 
support from lattice simulations and qualitative pictures. The reasons for quark confinement 
are somewhat complicated; no analytic proof exists that quantum chromodynamics should 
be confining, but intuitively, confinement is due to the force-carrying gluons having color 
charge. 

Asymptotic Freedom 



Asymptotic freedom which implies that the effective coupling constant = decreases 



logarithmically at short distances or high momentum transfers, a property which has a 
rigorous theoretical basis. This is the basis for perturbative QCD which is relevant for 
processes involving large momentum transfers, we have 

asm = .^^ ,\ (2.12) 

the running of as{Q'^) with Q'^.Aqcd is the QCD scale factor which effectively defines the 

energy scale at which the running coupling constant attains its maximum value. For |nj 

< 11, it is clear that asiQ"^) decreases as increases and approaches zero as Q"^ ^ oo or 

r — 7- 0. This is known as the asympotic freedom property of QCD |22j . 
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2.7.2 Quarkonia And Hadron Spectroscopy 

Quarkonium designates a flavorless meson whose constituents are a quark and its own an- 
tiquark i-e QQ, Q — c,b, is caUed quarkonium e.g. charmonium cc, bottomonium bb . 
Examples of quarkonia are the J/\E' and T{nS) . Because of the high mass of the top quark, 
a toponium does not exist, since the quark decays through the electroweak interaction before 
a bound state can form. 

For many reasons the strong interactions of hadrons containing heavy quarks are easier 
to understand than those of hadrons containing only light quarks. The flrst is asymptotic 
freedom, the fact that the effective couphng constant of QCD becomes weak in processes 
with large momentum transfer, corresponding to interactions at short distance scales. At 
large distances, on the other hand, the coupling becomes strong, leading to nonperturbative 
phenomena such as the confinement of quarks and gluons on a length scale Rhad ~ ^/ ^qcd ~ 
1/m, which determines the size of hadrons. Roughly speaking, Aqcd ~ 0.2GeV is the energy 
scale that separates the regions of large and small coupling constant. When the mass of a 
quark Q is much larger than this scale, ttiq ^ Aqcd, it is called a heavy quark. 

The light degrees of freedom arc blind to the flavor (mass). This is known as flavor 
symmetry. The heavy quark spin also decouples from the strong interaction. The decoupling 
of the spin in the heavy quark limit leads to the spin symmetry. These two symmetries have 
important consequences, especially for the decays of beauty hadrons to lighter hadrons. 
These symmetries are only true in the heavy quark limit and are violated at order ^"^^^ . 



Since quarks are fermions with spin 1/2, the wave function is antisymmetric with the 
exchange of particles Q and Q. Under particle exchange, we get with space coordinates 
exchange, a factor (—1)^, with spin coordinates exchange, a factor (—1)'^+^ and with charge 
exchange, a factor C (C is called C-parity). Hence Pauli principle gives 




(2.13) 
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In the spectroscopic notation, state is completely specified as 



n2^+% (2.14) 

where n is the principal quantum number and J is the total angular momentum. Thus for 
L = 0, we have the following states 

n C = +1, n = l,2,.... 
n^So C = -1, n = l,2,.... 

The ground state is therefore a hyperfine doublet 1^5*0(0^+) and 1'^ 5*1(1 ). For L = 1, we 
have the following states 

n^Pj J = +1 C = -l 
n^Pj J = 0,1,2 C = l 

Similarly we can write states for L = 2. 

It is noted that the state ^Di has the same quantum number as ^5*1, therefore they can 
mix. Most of these states have been discovered experimentally [22]. Some of the states are 
predicted, but have not been identified; others are unconfirmed. 

The computation of the properties of mesons in Quantum chromodynamics (QCD) is 
a fully non-perturbative one. As a result, the only general method available is a direct 
computation using Lattice QCD (LQCD) techniques. However, other techniques are effective 
for heavy quarkonia as well. The speed of the charm and the bottom quarks in their respective 
quarkonia is sufficiently smaller, so that relativistic effects in these states are much less. This 
technique is called non- relativistic QCD (NRQCD). 

All known hadrons are color singlets. The exchange of gluons can provide binding between 

quarks in a hadron. The two body one gluon exchange color electric potential is given by: 
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ks— ,ks = {-- for qq, 



(2.15) 




ks— ,ks = {-- for qq 



(2.16) 



Here i, j are flavor indices. Since the running coupling constant becomes smaller as we 
decrease the distance, the effective potential Vij in the lowest order as given by one-gluon 
exchange potential is a very good approximation for short distances. We conclude that 
for short distances, one can use the one gluon exchange potential, taking into account the 
running coupling constant ag- The second regime, i.e. for large r, QCD perturbation theory 
breaks down and we have the confinement of the quarks. One may look for the origin of 
this yet unsatisfactorily explained phenomena. There are many pictures which support the 
existence of a linear confining term. One of them is the string picture of hadrons. 

An early, but still effective, technique uses models of the effective potential to calculate 
masses of quarkonia states. In this technique, one uses the fact that the motion of the quarks 
that comprise the quarkonium state is nonrelativistic to assume that they move in a static 
potential, much like nonrelativistic models of the hydrogen atom. One of the most popular 
potential models is the so-called Cornell potential 



where r is the effective radius of the quarkonium state a and b are parameters. The first part, 
a/r corresponds to the potential induced by one-gluon exchange between the quark and its 
anti-quark, and is known as the Coulombic part of the potential. The second part, br the 
linear term is the phenomenological implementation of the confining force between quarks, 
and parameterizes the poorly-understood non-perturbative effects of QCD. Relativistic and 
other effects can be incorporated into this approach by adding extra terms to the potential. 



V{r) ^ - + br 



(2.17) 
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Chapter 3 



TETRAQUARKS: THE X (3872) 

3.1 Facts About The X (3872) 

The X was found in an exclusive decay in August 2003 [1]. 

B+ ^ K-^X (3872) ^ ir+ J/^vr+vr- (3.1) 

Belle measured its mass: 

m{X (3872)) = 3872.0 ± 0.6{stat.) ± 0.5 MeV/c^{syst.) (3.2) 
Belle also set a limit on its decay width: 

r(X (3872)) < 2.3 MeV^ (3.3) 

The most natural interpretation was a new coventional cc state. These exotic states refuse 
to obey conventional cc charmonium interpretation [21 13]. So many theories came: molecular 
models [41 [5], more general 4-quark interpretations including a diquark-antiquark model [6]- 
|10] . hybrid models [H] etc, to explain the nature of the X (3872) . In this chapter we mainly 
focus on the diquark-antidiquark model. 



3.1.1 Charmonium Hypothesis 

The fact that the X (3872) decays to J /ipTi^TT~ suggests that it contains c and c quarks, and 
the most natural choice to explain the X (3872) is an unseen charmonium state. Charmonium 
system has been thoroughly studied in [30]. After the discovery of the X (3872) more and 
more similar narrow resonances have been discovered and confirmed at electron-positron and 
proton-antiproton colliders. Twenty new unexpected charmonium-like particles have been 
found which have clear clashes with standard charmonium interpretations. Let us discuss the 
conventional charmonium states and the main objection for assigning them to the X(3872). 

The lightest charmonium state is called rjc- This is the S-state in which the spins of 
the quarks are antiparallel, so that the total spin is zero and the total angular momentum 
J = 0. The radial quantum number of the r/c is ?t- = 1, so that the spectroscopic notation {n 
"^^^^Ljpc) for this state is I^Sq and the J^^ is ^ . The next, a bit heavier, state is the J /ip. 
The next cc state Xco is a P— wave, more exactly I'^Pq) with J^^ = This is a part of a 
triplet, three particles with the same L and S, but different J. The other two particles in the 
triplet are Xci(l'^-Pi++) and Xc2(l^-f2++)- We discuss the upper part of the spectrum. Only a 
few of the charmonia states with masses above the DD threshold are cosidered. There are 
a few interesting states e.g. and ^D2. The decay of these states into DD is forbidden 
because of their spin-parity J^*" = 2~^. Both D and D have zero spin and the spin-parity 
of the DD system is determined by the same way as for two pions. The total even J of 
DD system constrains the C— and P— parities to be positive. Therefore these states cannot 
decay into DD and expected to have small widths. 

As mentioned earlier that initially the X (3872) was expected to be one of the so far 
unknown higher mass charmonium states. However, interpreting the X (3872) as a conven- 
tional state is problematic. We go through all the cc states which are not yet identified 
and evaluate by their suitability for the X (3872) . The states 25", 3P, 3D and higher are 
expected to be much heavy to associated with the X (3872) . We do not consider IS, IP 

and 25, because they are unambiguously identified already. Ten states remain, two of them 
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are known: l^Di — is '?/'(3770) and 3^Si — is '?/'(4040), so we will not consider them as serious 
candidates for the X (3872) . Now the eight states remain. The remaining eight states can 
not be interpreted as X [18], as mention in the following Table. To be precise, X (3872) is 
not a conventional charmonium state. 



Table 3.1: Conventional charmonium states and the main objection for assigning them to 
the X(3872). 



1 


2 


3 


4 




Mass 
MeV/c^ 


JPC 


Main objections for the X (3872) asignment 


13^2- 


~ 3838 
~ 3838 
~ 3838 


1 — 
0++ 
0++ 


expect rjcTin ^ J/tpmr 
not seen decay to Xci7 
not seen decay to Xc27 


23P2++ 


~ 3968 
~ 3932 
~ 4008 
~ 3966 


0++ 

1— 
1— 
1— 


wrong cos^^ j/^ distribution 
DD not suppressed— )-broad 
too large expected width to J/ip'j 
DD not suppressed— i-broad 


31^0-+ 


~ 4040 


1^- 


mass expected to be close to 



3.1.2 Weakly Bound D-D* State 

There are two possibilities to form bound states out of two quarks and two anti-quarks: 
In this configuration, binding each quark to an anti-quark [ga^"] and allowing interaction 
between the two color neutral pairs • Other possible configuration will be dis- 

cussed in the next section. The mainstream thought has been that of identifying most of 
these resonances as molecules of charm mesons. In particular the X (3872) happens to have a 
mass very close to the DD* threshold. The binding energy left for the X (3872) is consistent 
with E ~ 0.25± 0.40 MeV, thus making this state very large in size: order of ten times 
bigger than the typical range of strong interactions. These questions apply to other near-to- 
threshold hypothetical molecules and have induced to think to alternative explanations for 
the X (3872) and its relatives. 
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3.2 Diquark-Antidiquark Model 



In the diquark-antidiquark model the mass spectra are computed as in the non-relativistic 
constituent quark model. In the constituent quark model hadron masses are described by 
an effective Hamiltonian that takes as input the constituent quark masses and the spin-spin 
couplings between quarks. By extending this approach to diquark-antidiquark bound states 
it is possible to predict tetraquark mass spectra. The mass spectrum of tetraquarks QQ 
with Q = [qq] can be described in terms of the constituent diquark masses, uiq, spin-spin 
interactions inside the single diquark, spin-spin interaction between quark and antiquark 
belonging to two diquarks, spin-orbit, and purely orbital term [6]. This model has been 
tested on standard charmonia and has a rather good behavior to determine the mass spectra. 



3.3 Constituent Quarks and Spin-Spin Interactions 

In the costituent quark model the Hamiltonian is |50j : 



where the sum runs over the hadron constituents. The coefficient K-ij depends on the flavor of 
the constituents i, j and on the particular color state of the pair. Couplings for color singlet 
combinations are determined from the scalar and vector light mesons. For the L = mesons, 
taking us states, Eq. fl3.4p gives 




(3.4) 




(3.5) 



Mk = mg + m, + /C,g-[--] 
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Table 3.2: Constituent quark masses derived from the L = mesons and baryons. 



Constituent mass (MeV) 


Q 


s 


c 


b 


Mesons 


305 


490 


1670 


5008 


Baryons 


362 


546 


1721 


5050 



Table 3.3: Spin-Spin couplings for quark- antiquark pairs in in the color singlet state from 
the known mesons. 



Spin-spin couplings 


qq 


sq 


ss 


cq 


cs 


cc 




318 


200 


129 


71 


72 


59 



Similarly for the vector meson K* 

3 

Mk* = mg + m, + /C,g-[1(1 + 1) - 2] 

Mk* =mg + ms + ICggl^] 

Mk' -Mk^ 2]Csq, ICsq = 195 MeV 

Adding the similar equations for tt — p, D — D*, Ds — D*, J/ip — rjc complex we obtain the 
values of the spin-spin couplings, for quark-antiquark pairs in the color singlet state from 
known L = mesons. 

Now spin-spin coupling for quark-quark in color 3 (antitriplet) state can be calculated 
from the known L — baryons. The qq couplings are determined from the masses of the qqq 
baryons ground (J = 1/2) and excited (J = 3/2) states. Lets take the uds states:A, E, E*, 
which gives 

M = 2m, + m, + (/C,,), [3(3 + 1) - ^] + (/C,.) JJ( J + 1) - S{S + 1) - ^] (3.6) 
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Table 3.4: Spin-Spin couplings for quark-quark in color 3 state from the known baryons. 



Spin-spin couplings 


qq 


sq 


cq 


cs 




98 


65 


22 
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Writing equations for P, A+, involving only (/Cgg)^, and for Ac, Ec E* , involving (JCqg)^ and 
(JCgc)^- Also we consider the three Sc states which give (/Cgc)^ and (/Csc)^ couplings. These 
are given in Table 3.4 where one can see that the coupling strength decreases with increasing 
mass. 

It is observed that the diquark correlation decreases when one of the light quarks is 
strange. According to one gluon exchange (c.f Eq. fl2.15l) and Eq. fl2.16l) ). we have 

i^ij)3 = ^{^ij)o (3-7) 

This relates coupling of antitriplet to the singlet state. 

The couplings corresponding to the spin-spin interactions have been calculated for the 
color singlet and color antitriplet only. The couplings are not necessarily in the singlet 
state but octet couplings {ICcc)g are also possible. The quantities Kqq, fCcq and /Ccc involve 
both color singlet and color octet couplings between the quarks and antiquraks in a. QQ 
system (A quark in the diquark Q could have a color octet spin-spin interaction with an 
antiquark in the antidiquark Q). For the diquark attraction in the 3— color state, we can 
write = [eg]* = e^^^cjqk, where i, j, k are color indices in the fundamental representation 
of SU{3). The color singlet hadron is written as 

[cq][cq] = e'^^e,j,k'{cjqkWq^') = {cjC^){qkt) - {cjq^){qkC^) (3.8) 
Rearranging the color indices in the last term by using SU {N) identity for the Lie algebra 
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generators 



KjKi - 2(5i«5ifc - j^SijSki) 



(3.9) 



a=l 



where N is the number of colors. A color octet (A^ = 3) qq state can be written as q^X'^jq-' , 
and hence 



a=l 



i(lk){c' q ) - —{cjC>){qkq 



(3.10) 



Using Eg. ( 13. 9 P and Eq.f l3.10p . we extract the octet term as follows: 



[cq][cq] = {Cjc'){qkq)- 



1 1 



-[eX'l^c>){q^\l^q^) + -{cfc>){q,q 
1 



(3.11) 
(3.12) 



This formula gives information about the relative weights of a singlet and an octet color 
state in a diquark-antidiquark picture. We have three colors running in the sum qc* whereas 
a = 1, 8 in c^X^-cK Therefore the probability of finding a particular qq pair in color singlet, 
for example cc in the color singlet state CjC? , is half the probability of finding the same pair 
in color octet cA"c i-e, 3 X 2/3 = 1/2(8 x 1/2). 
We write for fCcc [6]: 

1 9 

(3.13) 



^c-c {[cq][cq]) = - (/Cec-)o + 2 (^cc-)8 



where (/Ccc)o is reported in Table 3.3. {}Ccc)s can be derived from the one gluon exchange 
model by using the relation [7] : 



(/Cee-)x - (X) - (3) - (3)) 



(3.14) 



where X is the color representation of the two quark system, with (X) = 0, 4/3, 4/3, 3 
for X = 0, 3, 3, 8 respectively. It is found that 
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(/Ccc)o ~ -g (^cc)8 ~ 3 = -g(^cc)o (3.15) 
Finally, from Eq. (l3.13p . one has 

^cc{[cq][cq]) = \ (^cc)o 

Now we have all the couplings and let us apply it to calculate the mass of light diquark for 
a simple case of ao(980): 

Co (980) = [sg]s=o[sg]s=o 
Using the Eg. (13. 5 p and eigen state given in the above equation, we calculate the mean value 



(ao ao) = Im^sq] + {^sc^A-^ ~ 2^ ^ - 3(/Cs<j)3 

m[,q] = 595 MeV 

Similarly one can calculate the m\^^^ = 396 MeV from (t(481). 

3.4 Spectrum Of Hidden Charm Diquark-antidiquark 
States 

In the diquark-antidiquark model effective Hamiltonian takes the form: 

H = 2mQ + Hfs^^ + Hfs^^ + HsL + HLL (3.16) 

where mg is the mass of diquark, -ff^f is the spin-spin interaction inside the single diquark, 

H^sT'^ is the spin-spin interaction between quark and antiquark belonging to two diquarks, 
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HsL is the spin-orbit, and is purely orbital term [6] i.e 



H'i's^ = 2(/C,,)3[(S,-S,) + (S,--S5)], (3.17) 

-^ss^^ = 2(/Ccg)(Sc ■ + Sc ■ Sg) + 

2/Ce5(S, ■ S,-) + 2]Cgg{Sg ■ Sg), (3.18) 

HsL = 2Aq{Sq-L + Sq-L), (3.19) 

H,L = Bs^^^^. (3.20) 



The overall factor of 2 is just a convention used in the literature. Aq, Bq are coefficients to 
be calculated by using known data. We will use these values in the next chapter. 

To calculate the spin-spin interaction of the QQ states explicitly, we use the following 
non-relativistic notation for labelling the state 

\Sq, Sq, J) = |r, r; J) = {c^Ta,u''){cT'^,u') (3.21) 



where, Sq and 5*2 are the spin of diquark and antidiquark, respectively, J is the total angular 
momentum and the are 2x2 matrices in spinor space. Using Pauli matrices these can 
be written as: 

for spin and 1, respectively. The matrices F are normalised so that: 

Tr[(F'^)t(F^)] = 



We define the spinor operators as: 



SJF) 



Scir) 



-a 



^F 
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since, 

T 

We calculated formula for total spin operator as expected: 



(S„+Se)|rO> = (3.23) 

[(S J+( SJ^] |F> = ie^^^\Y^) (3.24) 

We also find: 

(0|S„|1) = -(0|Se|l) = ^; (3.25) 

(1|SJ1) - (1| Se|l) = ^(l|(S„+Se)|l) (3.26) 

We have used the following Pauli matrices properties: 

Tr(a'a') = Tr(I) ^ 2 

By using this information we now calculate the matrix elements of products of spin operators. 
There are two cases. 

Same diqueirk, e.g.Sy, ■ Sc- 

This operator is only a combination of Casimir operators and is diagonal in the basis. 

2(S„-S,) = (Se„)'-(SJ^-(SJ^ (3.27) 
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Different diquarks, e.g.Su ■ 5"^^. First we consider J = states, represented by 



|0s,02; Oj) 

|1q,1q; Oj) 

Using the basic definitions, we get: 

2(S„- S^)|OQ,Og; Oj) = 
2(Su ■ Su) |1q, 1q; Oj) = 



- ((72) ® (as) 



2^3 



(0-20-*) ® (0-20-*) 



1 /3 

- ((T2(T^) ® (0-2(1*) = ^ |1q, Ig; Oj) (3.28) 

Oq,Oq; Oj)-|1q,1q; Oj) (3.29) 



4^3 
2 



wliicli leads to tlie following matrices: 



2(5*11 ■ Su] 



^ 2 



\ 2 



Now we consider J = 1 states, given in the tensor basis: 



(3.30) 



|0q,1q; 1j) 
|1q,0q; Ij) 



^ ((12) ® ((^2(7') 
^ ((72(7*) ® ((T2) 



2^2 



r'^' ((T2(T^) ® (a2(x') 



(3.31) 



The normalisation of the Hamiltonian in Eg. (13. 161) . using the basis of states defined above 
with definite diquark and antidiquark spin and total angular momentum, will give the spec- 
trum of diquark-antiquark states. There are two different possibilities: Lowest lying [cq] [cq] 
states {Lqq = 0) and higher mass [cq] [cq] states {Lqq = l) . 
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3.4.1 Lowest Lying [eg] [eg] States 

In the ground state the two diquarks interact only by spin couphngs because the angular 
momentum is zero {Lqq — 0). An effective non-relativistic Hamiltonian can be written 
including spin-spin interactions within a diquark and between quarks in different diquarks. 
The states can be classified in terms of the diquark and antidiquark spin, Sq and Sq, total 
angular momentum J, parity, P and charge conjugation, C. Considering both good {Sq — 0) 
and bad {Sq = 1) diquraks and having Lqq = we have six possible states which are listed 
below. 

i. Two states with J^^ = 0++: 



|0++> = |OQ,Og; Oj); (3.32) 

|0++'> = |ls,lg; Oj). (3.33) 

ii. Three states with J = 1: 

|i++> = -^{IOqAq; 1j) + \1q,0q; ij)); (3.34) 

|l+-> = ^(|0q,1s; 1j)-|1q,0s; U); (3.35) 

|l+-'> = \1q,1q, Ij). (3.36) 



All these states have positive parity as both the good and bad diquarks have positive parity 
and Lqq — 0. The difference is in the charge conjugation quantum number, the state [l"'"'') 
is even under charge conjugation, whereas ll"*"") and ll"*""') are odd. 
iii. One state with J^*^ = 2++: 

|2++> = |lQ,lg; 2j). (3.37) 

Keeping in mind that for Lqq = there is no spin-orbit and purely orbital term, the 
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Hamiltonian of Eq.f l3.16|) takes the form 



H = 2mQ + H?,^^ + Hif^ 



(3.38) 



Thus, 



H 



2m[g,] + 2(/C,,)3[(Sc ■ S,) + (S,- ■ S,-)] + 2/C,,-(S, ■ S,-) 



+2(/Ccg)(Sc ■ Sq + Sg ■ Sg) + 2/Ccc(Sc • Sg). 



(3.39) 



The diagonahsation of the this Hamiltonian with the states defined above gives the eigen- 
values which are needed to estimate the masses of these states. It is calculated that for the 
l"*""*" and 2++ states the Hamiltonian is diagonal with the eigenvalues 



All other quantities are now specified except the mass of the constituent diquark. The 
X(3872) is a [eg] 5=1 [eg] 5=0 tetraquark. By diagonalizing the Hamiltonian in Eq.f l3.39p and, 
using the spin couplings derived above, the mass of the diquark [eg] was fixed by using the 
mass of X(3872) as input, yielding m[cg] = 1.933 GeV. The 1^^ state is a good candidate to 
explain the properties of X(3872). In order to reduce the experimental information needed 
we estimate the remaining diquark masses by substituting the costituent quark forming the 
diquark. We have 



M (1++) = 2m[cq] - iJCcqh + -j^qq " ^cq + -^K^cc, 
M(2++) = 2m[e,] + (/Ce,)3 + ^/C,,- + /Ce,- + ^/Cec-. 



(3.40) 



(3.41) 



^[cq] -mq + m, 



'S 



(3.42) 



m[bg] 



m[cq] -rric + mh 



(3.43) 



mibg] -nig + m, 



'S 



(3.44) 



36 



Table 3.5: Diquark masses 



Diquark masses 


[ud] [sq] [cq] 


[cs] [hq] 


MeV 


396 595 1933 2118 


5119 



Now, we have all the input parameters to calculate the mass spectrum numerically. 
Putting masses of diquarks from Table (2.4) and values of couplings from Tables (2.2), (2.3) 
in Eq.f l3.40l) . we get the mass for the hidden cc tetraquark 1++ state: 



M (1++) = 3.872 GeV 



(3.45) 



Taking the X (3872) as input we can also predict the existence of a 2^^ state that can be 
associated to the X (3940) observed by Belle 



M (2++) = 3.952 GeV. 



For the corresponding 0"'""'" and 1^~ (Labelled as Z ) tetraquark states, the Hamiltonian 
is not diagonal and we calculated the following 2x2 matrices, using the non-relativistic 
notation for labelling the state as in Eq. fl3.30p : 



M (0++) 



— 3(A^cg)3 ^ i^qq + ^cc ~ '^^cq) 



\ 



^ ^ i^qg ~^ ^cc ~ 2A^cg) (^05)3 ~" i^qq + ^cc + '^^cq) y 



M(l+-) 



v 



-(M3 + /C,,-- 



^qq 



^qq ^cc 

(IT \_ IT (K.qq+K.cc.) 

yl^cqjZ l^cq 2 



To estimate the masses of these two states, one has to diagonalise the above matrices. After 
doing this, the mass spectrum of these states is shown in Fig 3-1. 
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Figure 3.1: Lowest lying hidden charm spectrum. 



3.5 Isospin Breaking and Decay Widths of the J^^ = 



1++ Tetraquark 



3.5.1 Isospin Breaking 

In this section we discuss the isospin breaking effects which were neglected in the previous 
section. The isospin quantum number is related to the finer structure of the X state. The 
two flavor eigenstates and Xd mix through self energy diagrams, which annihilate a 
uu pair and convert it into a (id pair through intermediate gluons. In the basis {X„,Xrf} 
the annihilation diagrams contribute equally to all the entries of the mass matrix, while the 
contribution of the quark masses is diagonal. The resulting 2x2 mixing mass matrix is: 



V 



2mu + 5 5 
5 2md + 6 



J 



(3.46) 



where S is the contribution from quark annihilation diagrams. At the scale determined by 
the cc pair the annihilation term 6 is expected to be small and thus the mass eigenstates 
should coincide with flavor eigenstates to a rather good extent. Isosp in-breaking introduces 
a mass splitting and the mass eigenstates called X[c,i] and X[c,h] (for lighter and heavier of 
the two) become linear combinations of X^^cu] and X^cd]- One can put: 
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X[c,i] = COS 6 X[cu] + sin 6 X[cd] 
X[c,h] = - sin 6* X[cu] + cos 6 X[cd] 



(3.47) 
(3.48) 



The mass differences are estimated to be small, where ^ is a mixing angle. The electromag 
netic couplings of the tetraquarks X^cfl and X^^ch] will depend on the mixing angle 9. 



We get, 



X[c,l] ^ 




X[cu] ' 













COS 6 sin 6 
— sin 9 cos 9 



2m„ = M{X[,,i])cos^9 + M{X[,,h])sm^9 
2md = M(X[,,;])sin2^ + M(X[,,^])cos2^ 



Take the difference: 



M(X[,,^]) - M(X[e,,]) = 2(m, - m„)/ cos(2^^) (3.49) 

Infact, there are two different states X(3872) and X(3875) which were not excluded from the 

experimental data [ITJ HH]. The isospin violation in the tetraquark picture is the possibility of 

u — p^ mixing, as proposed in [D]. The observation in 2006 of a state decaying to D^D^n with 

mass 3875 MeV favored the assignment: Xu = X(3875), decaying mainly into J/ipn'^Tc' and 

Xd = X(3872) decaying into D^D^n. The mass ordering of these two neutral states seems 

to be reversed, since the u quark is lighter than the d quark and thus one would expect X„ 

to be lighter than X^. However the quarks which form the diquarks in the X^ have the same 

electric charge and thus a consistent consideration of the electrostatic energy can perhaps 

change the order of the masses. Besides these two neutral states, two charged states arise as 
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a natural prediction of the tetraquark picture X+ = [cu] [cd\ and X = [cd] [cu] . The charged 
partners are X"^ are not observed [3T] . 

3.5.2 Decay Widths of the J^^ = 1++ 

Originally the X(3872) was found through its decay into J/iljTr~^TT~ , but other decay modes 
were also investigate. The decay of a diquark-antidiquark bound state into a pair of mesons 
can occur through the exchange of a quark and an antiquark belonging respectively to the 
diquark and the antidiquark. There are indeed three different flavor configurations. Thus 
we need to introduce three amplitudes. Two of them account for X — )■ D^D^*: 

Third one the exchange of a light quark and a heavy quark accounts for, the charmonium 
channels i-e, 

A{[cq][cq] ^ [qq][cc]) = A (3.50) 

The only available ones are J/ip ^ 2tt and J/ip Svr, dominated by and u respectively 
and were confirmed experimentaly. 




fcj X — ^ yy 

Figure 3.2: Radiative and Hadronic decay of the X(3872) described with the same contact 
vertex A. The radiative decay proceeds through the hadronic transition. 



3.5.3 Hadronic Decays 

The decay rate for X — J/?/' + /, as shown in Figure 3-2 (a) can be written as 
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ds 



SttM 



MvT 



TT {S- M2)2 + [MvTvy 



(3.51) 



with 



p the decay momentum: 



/ = TT^TT (tt^tt Tr*^) for V = p{u) 



pis) 



^A(Mx[,,j, M^, Mv) 



Id] 

-2{M^Mvf 



(3.52) 



(3.53) 



where the coefficient xi y is: 



{cos 6 ±sin^)" 



and A is taken to be 2.6 GeV [9\. Similarly we can derive all above equations for higher 
mass state Xic,h]- Numerical integration of 



n 



p{s) 



[s - M2)2 + (M,r, 



ds 



p{s) 



(3?Ti^)2 



[s - M2)2 + (M^r, 



gives 



{p)p = 126 MeV, {p)^ = 22 MeV 
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r(X[e,i] ^ ^ + TT+TT-) 



id] 



r(X[,,^]^V' + 7r+7r-7r°) = 



0, _ 2x;,^ |AP 



c,(] 



r(X[c,i] ^ -0 + TT+TT") = 2xi,p • 2.3 Mey = 3.78 MeV (3.54) 

r(X[c,i] ^ + tt+tt-tt") = 2xi,^ ■ OAMeV = 0.66 MeV (3.55) 

Using these values of decay rates, we can get the information about the mixing angle 9. 



M = (^2il±ii^i^.(% = 0.802 (3.56) 
(EfM) - (COS 9 -sinO)^ W-^-0S02 (357) 

where from the Belle experiment we have 

Vr(27r)^,,, 

Putting every thing together we have 9 = ±20° , for X^cfi and -^[c,/i] respectively. For the 
charged state X^, that decay via yO-exchange only: we have 



3.5.4 Radiative Decays 

The amplitude for the radiative decay X ^ J/t/j^ proceeds through the annihilation of a 

pair of light quarks into a photon, the hadronic part of the amplitude is the same as in 
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the decay X — )■ J/iIjtt~^tt~ . The radiative decay proceeds through the hadronic transition 
X — )■ J/ 1pp. Exploiting the Vector Meson Dominance (VMD), which describe interactions 
between photons and hadronic matter [32], one can write transition matrix from the Fig. 
3-2 (b) as: 



(J/^7|X) = {l\p)—,{J/MX) = (3.58) 
nip nip 



Thus the partial decay width is: 



r(.Y ^ j/in) = 2iA'i ( I^) ' ' VMMv. M,. 0) ^^^^^ 



P 



ml) 8nMj. 2Mx 



Using fp = 0.152 GeV^ ymj we get, 

T{X J/t/j'y) 



T{X J/V'vr+TT-) 

which is in agreement with experimental value [28] . 

T{X ^ J/ip-f) 
r{X J/ijjTr+n-) 



0.44 



< 0.40 



Similarly it is easy to calculate other radiative decays. The experimental values for these 
radiative decays are [28] : 



r(X ^ J/ljjTT+n-) 

TjX ^ Xc27) 
r(X ^ J/-iPn+n-) 



1.1 



We exploit the result obtained for the width of X — )■ J/ip'j to give an estimate of the 
decay width into J/ip'j'j. We compute the transition matrix element in terms of A using the 
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coupling of the J/ip to the 7 for the Figure 3-2 (c). 



(77|X) = (7| J/^)^(J/^7|X) = ^(J/^7|X) = M^^A (3.60) 



Thus the partial decay width is: 



r(X^77)4|A^I^ f4V^^^^^^^ (3.61^ 



3 \^J/^/ V'^p/ 8^^i 2Mx 



Using fj/^ = 1.254 GeV^ [ID] we get, 

r(x -> 77) 



3 X 10 



-5 



r(X ^ J/^TT+vr-) 
that is greater than the upper limit provided in experimental data 

r(X^77) ....,..-5 



r(x ^ j/^Ti+Ti-) 



< 1.5 X 10" 



The inconsistency of the theoretical prediction with respect to data is not dramatic if we 
take into account the very strong assumptions made to derive the Eq. (l3.6ip . 
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Chapter 4 

GROWING EVIDENCE OF 
TETRAQUARKS: THE ^^(lOSQO) 

4.1 Introduction 

Using the diquark-antidiquark model, in the previous chapter we learnt to derive the spec- 
trum for lowest lying [eg] [eg] states with q = u, d, and discuss the decays of X(3872). 
This model can be easily applied to lowest lying [6g][&g] states with q = u, d, s, and e. 
In this chapter we will discuss the spectrum of the higher mass [&g][&g] tetraquark states 
with q = u, d. Particularly we are interested in the decays of higher tetraquark particles 
with q = u, d. First, there is evidence for ss bound state, Ys{2n5) having the quantum 
numbers J^*" = 1 , first observed by BaBar in the initial state radiation [ISR) process 
e"'"e~ — )■ 7isR /o(980), 0(1020), where /o(980) is the O"*""*" scalar state [33]. This was later 
confirmed by BES [34] and Belle [35]. In December 2007, the Belle collaboration working at 
the KEKB e+e~ collider in Tsukuba, Japan, reported the first observation of the processes 
e"'"e~ — )■ Yibq] — )■ T(1S', 25") vr+vr" near the peak of the T(5S') resonance at the center-of-mass 
energy of about 10.87 GeV [12j. Belle measurements near the T(5S'), however, did not fall 
in line with theoretical expectations [12]. Their data were enigmatic in the partial decay 
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widths for T{5S) — )■ T(IS') tt+tt^ and T{2S) tt+tt^ which were typically three orders of 
magnitude larger than anticipated in QCD [36]. Production and decays of the T{nS) states 
(n being the principal quantum number) are popular theoretical laboratories to test QCD. 
In particular, the final states T(1S', 25*) tt+tt^ arising from the production and decays of the 
lower bottomonia states, such as T(45') — > T(15') tt+tt^, have been studied in a number of 
experiments over the last thirty years and are theoretically well-understood in QCD 
In addition, the dipion invariant mass distributions in these events were distinctly different 
from theoretical expectations as well as from the corresponding measurements at the T(4S'), 
undertaken previously by Belle. 

In the conventional Quarkonium theory, there is no place for such a nearby additional 
bb resonance having the quantum numbers of T{5S). An important issue is whether the 
puzzling events seen by Belle stem from the decays of the T{5S), or from another particle Y/, 
having a mass close enough to the mass of the T{5S). We will see that the interpretation of 
the Belle data is that the anomalous T(1S', 25*) tt+tt^ events are not due to the production 
and decays of the T{5S), but rather from the production of a completely different hadron 
species, tetraquark hadrons with the quark structure Y^^g] = [bq]\bq] states with q = u, d 
and their subsequent decays. A. Ali et. al. call this state a "Brand New Form of Matter". 
Identifying the J^^ = 1~~ state (10900) seen in the energy scan of the e~^e bb cross 
section by BaBar |31] with the state y[6g] (10890) seen by Belle [12]. 

Clearly, two aspects of the Belle data had to be explained: 

(a) the anomalously large partial decay rates and 

(b) the invariant mass distributions of the dipions. 

A dynamical model based on the tetraquark interpretation of Yj^g] (10890) was presented in 
[38] where it was pointed that it is in agreement with the measured distributions in the decays 
Y[bg] — 7- T(IS') TT+TT^, T(25') TT+TT^. They have argued that the decays Y^^g] — t- T(15', 28)71^11' 
are radically different than the similar dipion transitions measured in the T{AS) and lower 
mass Quarkonia. Most resently the decays Y^hg] — )■ T{1S)K^K' was investigated by A. Ali 
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et al. [39] further supporting the Belle data. We will discuss in detail the anomalously large 
partial decay rates and not the invariant mass distributions of the dipions. 



4.2 The Fifth Quark Flavor: Bottom Mesons 

Fifth quark was discovered, when in 1977 the upsilon meson ) was found 

experimentally as a narrow resonance at Fermi Lab. with mass ~ 9.5 GeV. This was later 
confirmed in e^e~ experiments at DESY and CESR which determined its mass to be 9460±10 
MeV and also its width. The updated parameters of this resonance are mass 9460.37 ± 0.21 
MeV and width 52.5± 1.8 keV. Again the narrow width in spite of large phase space available 
suggests the existence of a fifth quark flavor called beauty, with a new quantum number 
B = —1 for the bottom (b) quark. With this assignment the formula Q = I3 + 1/2{Y + B + C) 
would give the charge of b quark the value — 1/3(/3 = 0). The mass of b quark is expected 
to be around 4.9 GeV as suggested by the T mass which is regarded as a ^Si bound state of 
bb. 

One would also expect the particles with B = ±1, such as bq or qb. The lowest lying 
bound states bq and qb have been found experimentally. The B = —1 states {B^, B~)Bg 
form an SU(3) triplet (3) and B = +1 states {B^, B^)B^ form another triplet (3). For 
p-wave multiplets 



{qb)L=i J^ = 2+,r 



JP = l+,0^ 



f D*0 75*0 



i=3/2 



^ J, = 1/2 

The masses and decay time of B-mesons are given below 



B^ = 5279.16 ±0.31 Mel/, r = (1.638 ± 0.11) x 10"^^ sec 
5° = 5279.53 ± 0.33 MeV, r = (1.530 ± 0.069) x lO"^^ 
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4.3 Spectrum Of Higher Mass [6g][6(/] States 

For the orbital excitation Lqq = 1 we consider both good and bad diquarks. The orbital 
excitation Lqq = 1 leads to negative parity states, where 1 multiplet is of main interest in 
this chapter. To estimate the masses, we repeat the diagonalization with the basis: 



|1) = |0q,02; Ij) 

|2) = ^(|0q,1s; 1j) + |1q,02; Ij)) 

|3) = Ij) (4.1) 



Since for the both good and bad diquarks parity is positive as from Eq. fl2.8p and Eq. fl2.9p 
(0+ and 1+ respectively), the state |2) has P = C = —1, provided that Lqq = 1. Since 



therefore for the states |1) and |3) , Cqq = — 1 provided that Sqq = 0, 2 and Lqq = 1. First 
of all we have to change the Hamiltonian given in Eq. (13.161) by replacing the charm quark 
by bottom quark. 

H = 2mQ + Hfs^^ + Hfs^^ + Hsl + Hll (4.2) 

where: 

= 2(/C,,)3[(S,-S,) + (S5-S,)], (4.3) 
^ 2(/Cfe,-)(S6-S5 + S5-S,) + 2/C,5(S,-S5)+2/C,,(S,-S5), (4.4) 
Hsl = 2Aq{Sq-L + Sq-L), (4.5) 
Hll - B,^. (4.6) 
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To perform the digonalization we use the same shorthand notation as described in the previ- 
ous chapter (c.f. Eg. (13.2 ip to Eg. (13.301) ) for the basis vectors defined in Eg. (l4.ip . We derive 
the mass term shift Amss for higher mass [bq] [bq] states, due to the part of the Hamiltonian 
containing only spin-spin interaction terms, Hss- Let us first consider 



55 J 33 



Here we use |[sq(sq + 1) — |] = Sq.Sq 



(Hss) 



33 = 0-qAq', lj|-f^55^^|ls> 1q; 1j) + (1s>1s; '^j\Hs's'^'\Iq,Iq', Ij) (4.7) 
= 2 (/C,,)3 (2 - 3/2) + 2/C,,-[(2 - 3/2)^ + (2 - 3/2)^] 

+2/C,,-(2 - 3/2)1 + 2/C,5(2 - 3/2)i 

= (^69)3 - - -ICgg - -JChb (4.8) 



(SS)| 



Similarly we can easily calculate (Hss) 11 and (Hss) 22- AH off diagonal elements are zero. 
Thus finally we have. 



Am 



55 



-3 {}Cbq)l 












+ {JCqq + /C,5) /2 (4.9) 
Q Q i^bq)3, - ^bq 

\ _(/C^_ + /C,5)/2 I 

The eigenvalues of the spin-orbit and angular momentum operators given in Eg. (14. 5 p and 
Eg. (14.61) . were calculated by Polosa et al. we have 
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Aq = h MeV, for q = u, d, 
Bq = 408 MeV, for q = u, d, 



We use these values in order to calculate the numerical values of these states. Hence the 
eight tetraquark states {q = u, d) having the quantum numbers 1 are: 



y\bq] 


(Sq 


= 0, 


Sq 


= 0, Sqq 


= 0, ^22 = 1) 


= 2m[bg] + Ai + Bq, 


(4.10) 


y\bq] 


(Sq 


= 1, 


Sq 


= 0, Sqq 


= 1, Lqq = 1) 


= 2m[bg] + A + A2 - 


2Aq + Bq, (4.11) 


Y[bq] 


{Sq 


= 1, 


Sq 


= 1' Sqq 


= 0, Lqq = 1) 


= 2mibg] + 2A + A3 - 


4-^2, (4.12) 


Y[bq] 


{Sq 


= 1, 


Sq 


= 1; Sqq 


= 2, Lqq = l) 


= 2m[bg] + 2A + A3 - 


-GAq + BqIAAS) 



where A's are the diagonal elements of the matrix AMss given in Eq. fl4.9p . The quantity A, 
is the mass difference of the good and the bad diquarks i.e. 

A = mQ {Sq = l)-mQ {Sq = 0) . 

and from Eq. (13.431) 

m^bg] = m[cq] -mc + mb 

Now one of the remaining unknowns in this calculation is the quantity A, the mass difference 
of the good and the bad diquarks. Following Jaffe and Wilczek [I9|, the value of A for diquark 
[bq] is 202 MeV. 

We recall previous chapter where we have used the known mesons and baryons to calculate 
the couplings of the spin-spin interaction. We can extend the same procedure to the 5 = 1, 
L = (0, 1) meson states B*, Bi (5721), B2 (5747) to calculate the values: 

Ultimately by putting things together, the masses for the states given in Eqs. (14. 10114.131) 
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Table 4.1: Spin-Spin couplings for quark-antiquark pairs in in the color singlet state from 
the known bottom mesons. 



Spin-spin couplings 


qq 


hq 


hh 


(/C,,)o(MeV) 


318 


23 


36 



Table 4.2: Spin-Spin couplings for quark-quark in color 3 state from the known bottom 
baryons. 



Spin-Spin couplings 


qq 


hq 


(/C.,)3(MeV) 


98 


6 



are given in Table and can be compared with the ones estimated in refs. ^0] using the QCD 
sum rules. 

Note that there are 8 electrically neutral self-conjugate 1 tetraquark states Y^^^j'* with 
the quark contents with q = d oi which the two corresponding to [6m] [few] and 

i.e., yj["j and Yj^^^l are degenerate in mass due to the isospin symmetry. Their 
mass difference is induced by isospin splitting — mu-, mixing angle 9 and is estimated as 
AM(Y5) = (5.6±2.8) MeV. Due to this small differnce in the following we will not distinguish 
between the lighter and the heavier of these states and denote them by the common symbol 
Yi). There are yet more electrically neutral J^^ = 1 states with the mixed light quark 
content [6c?] [6s] and their charge conjugates [6s] [M]. However, these mixed states don't 
couple directly to the photons, or the gluon, and are not of our main interest. 

Table 4.3: Masses of the the 1 neutral tetraquark states My"'* in GeV. The value My(i) 

l*"?! [bq] 

(for q = u,d) is fixed to be 10.890 GeV, identifying this with the mass of the Yh from BELLE. 







yibq] 






q = u,d 


10890 


11130 


11257 


11227 
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4.4 Decay Widths Of The ^6(10890) 



4.4.1 Leptonic Decay Widths 

For bottomonium systems, the corresponding decay widths are determined by the wave 
functions at the origin for the T{nS), \I/ft5(0), and by the derivative of these functions at the 
origin, ^''^^(O), for the P-waves. To take into account the possibly larger hadronic size of the 
tetraquarks compared to that of the bb mesons, we modify the Quarkonia potential, usually 
taken as a sum of linear (confining) and Coulombic (short-distance) parts. For example, the 
Buchmiiller-Tye QQ potential has the asymptotic forms 

Vir) ~ kggr , (for r -> oo) (4.14) 

V{r) ~ l/rln(l/A|cD , (for r ^ 0) (4.15) 

where kqQ is the string tension and Aqcd is the QCD scale parameter. The bound state 
tetraquark potential Vgg(r) will differ from the Quarkonia potential Ygqlr) in the linear part, 
as the string tension in a diquark fcgg is expected to be different than the corresponding 
string tension kqq in the QQ mesons. The diquarks-antidiquarks in the tetraquarks and 
the quarks-antiquarks in the mesons are in the same 3c3c color configuration, the Coulomb 
(short-distance) parts of the potentials will be similar. Defining, 

K=kQQ/kQQ (4.16) 

we expect k, to have a value in the range k G [|, ^] [42J. This will modify the tetraquark wave 

functions "^QQiO) from the corresponding wave functions of the bound bb systems, effecting 

the decay amplitudes and hence all the decay widths of the tetraquarks. The corresponding 

value for the tetraquark states [bq] \bq] is then calculated taking into account the ratio of the 

string tensions k. As the linear part of the confining potential determines essentially the 
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heavy Quarkonia wave functions, we find that to a good approximation: 



^ss(O) ^ ^^tm (4.17) 

which is what we can use in our derivations of the decay widths. 

The partial electronic widths ree(^[f,,i]) and Tee(Y[ti^h]) are given by the well known Van 
Royen-Weisskopf formula for the P-states, which we write as: 



167rg2a2|^'^ 12 



2 



(4.18) 



where Q = —2/3 is the diquark charge in Ytd = [bd][bd] and Q = +1/3 is the charge of the 
diquarks in Ybu = [bu]\bu], a = 1/137 is the electromagnetic couphng constant to lowest order 
and ^gg(^ = 'ip{4', 0)R'{r) is the first derivative in r of the wave function of the tetraquark, 
which needs to be taken at the origin i.e., 



*Qq(-1 = \/|^^'(0) (4-19) 



The wave function ^l/'j,^(0) depends on the underlying potential model. Tetraquarks are bound 
state of diquarks-antidiquarks with the same internal color structure as the Quarkonia, i.e. 
both are (3, 3) bound states. Hence, they also depend, in principle, on the potential models. 
We will not distinguish between the lighter and the heavier of these states because their 
mass difference is induced by isospin splitting — m„ and a mixing angle and is estimated 
as AM{Yh) = (5.6 ± 2.8) MeV. So we denote them by the common symbol Yb. 

To calculate the radial wave function, we took the value calculated by using the QQ- 
onia package [27] yielding |i?'(0)r = 2M2GeV. Since each derivative increases the energy 
dimension by and thus we need to normalize by kinetic energy u ^ rriQ of the diquark. For 
the lowest lying 1 state we get. 
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Tee ~ 0.12KeV (4.20) 

Since all the 1 states are P- waves, the -R'(O) value will not change because the masses of 
the diquarks remain the same. So the value of Fee only varies with the mass and therefore 
does not change significantly. This value is close to the experimental value given in 



4.4.2 Hadronic Decays 

We discuss the two-body hadronic decays of the yfe(10890), ^(g) Bi*\k)Bi*\l) . The 
decays Yj, — )• T(1S', 25") tt+tt^ are also Zweig allowed. These decays y[{,g] T{1S,2S) tt+tt" 
are too much phase space suppressed and require a dynamical model, which we will also 
discuss later. 

The Vertices of the dominant two-body hadronic decays of the y;,(10890) are [43j: 
Yb^ BB = Fik^" - P) 

Yh^BB* = f,e^'''"'kpl„ (4.21) 

Yh^B*B* = Fig^'Piq + iy -gf^'ik + qY + gf^iq + kY) 
The corresponding decay widths are respectively: 

^ r = 3S (4-22) 



r = 



F^|fc|3(48|A:|^-104M^|A:p+27M*) 



27r(M3-4|fc|2M)2 

The decaying momentum \k\ is as given in Eq. ( 13.521) 



_ y/M^ - (M, + - (Mfc - MiY 

|fc| - 1^ , (4.23) 

54 



where M is the mass of the decaying particle and M^, Mi are the masses of the decay 
products. The matrix elements are obtained by multiplying the vertices by the polarization 
vectors e^*^^. The polarisation vectors e^*^'^ satisfy the transversality condition with the 
polarisation sum 

= -g^'' + ^ (4.24) 

i ^ 

Thus, for the decay Yj^g] — )■ BqBq, the Lorentz- invariant matrix element is: 

M = eJ''lF(A;^ - P) (4.25) 

The decay constants F are non-perturbative quantities, which are beyond the scope in our 
approximation. We estimate them using the known two-body decays of T(5S'), which are 
described by the same vertices as given Eq.f l4.21|) [T6] . We use the decay widths for the 
decays T(5S') — )■ BB,BB*,B*B* from the PDG, yielding the coupling constants, called 
FpDG and \k\. 

Table 4.4: 2-body decays T{5S) B^*^B^*\ which we use as a reference, with the mass and 
the decay widths taken from PDG. The extracted values of the coupling constants F and 
the centre of mass momentum \k\ are also shown. 



Process 


TpDG [MeV] 


FpDG 


\k\ [GeV] 


T{5S) BB 


< 13.2 


< 2.15 


1.3 


T{5S) BB* 


15.4tn 


q 7+U.7 
' -0.9 


1.2 


T{5S) B*B* 


48tli 


1 +U.i3 


1.0 



The different hadronic sizes of the bb Onia states and the tetraquarks Y^^g] are taken into 
account by the quantity k, discussed earlier. 

Thus, the decay widths of Yj^Jj are consistent with the corresponding measurements 
by Belle. The other three higher 1 states have much larger decay widths and will be 
correspondingly more difficult to find. 
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Table 4.5: Reduced partial decay widths for the tetraquarks YjJ^J], the extracted value of 

the coupling constant F and the centre of mass momentum \k\. The errors in the entries 
correspond to the errors in the decay widths in Table 4-5. 



Decay Mode 


T/K^[MeV] 


F 


\k\[GeV] 




< 15 


2.15 


1.3 




18tl 


3.7 


1.2 


^ B*B* 


56t\t 


1 


1.1 



4.4.3 Dynamical Model For Yb ^ T(1S', 2S) tt+tt 
(Introduction and results) 

Explaining the larger decay rates for the transitions Yb{q) — )■ T{p) + 7c~^{ki) + 7T~{k2) was 
not so difficult, as the decays of Yb involve a recombination of the initial four quarks, as 
exemplified below by the process y[5„] = — ?■ {bb){uu), with the subsequent projection 
(bb) —7- T(15') and {uu) — )■ Tr+vr^.Such quark recombination processes do not require the 
emission and absorption of gluons, and are appropriately called Zweig-allowed, after the co- 
discoverer of the quark-model, George Zweig. The relevant diagrams for the decays Yi,{q) — )■ 
T(p) + 7r+(A;i) + 7r-(fc2). 

The decays Y^ — t- T(IS') tt+tt^, T(2S') tt+tt^ are sub-dominant, but Zweig-allowed and 
involve essentially the quark rearrangements shown below. With the J^*" of the Yh and 
T{nS) both 1 , the tt+tt" states in the decays Yh — t- T(IS') tt+tt", T(25') tt+tt" are allowed 
to have the 0++ and 2++ quantum numbers. There are only three low-lying states in the 
Particle Data Group (PDG) which can contribute as intermediate states, namely the two 
0++ states, /o(600) and /o(980), which, following [lH |15] we take as the lowest tetraquark 
states, and the 2++ gg-meson state /2(1270), all of which decay dominantly into tttt. For the 
decay Yf, — )■ T(15') tt+tt", all three states contribute. However, kinematics allows only the 

/o(600) in the decay Yb T(2S') tt+tt". This model encodes all these features. Finally the 
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fits of the Belle data with this model, yield the measured partial decay widths Ty{is)+2tt 
0.59 ± 0.04 ± 0.09 MeV and rT(25)+2^ = 0.85 ± 0.07 ± 0.16 MeV [38]. 



4.4.4 Radiative Decays 

The amplitude for the radiative decay Yf, — >■ xtj will proceed through the annihilation of a 
pair of light quarks into a photon. 




Figure 4.1: Radiative decay of the y(10890). The radiative decay proceeds through the 
hadronic transition. 

Exploiting the Vector Meson Dominance (VMD), one can write transition matrix from 
the Fie. 4-1 as: 



Thus the partial decay width is: 



2, I fpY 1 VHMn, 0) 



p/ Yt J 6 

where ^'^'^^^^My^^^''' ""^ decay momentum and the value of A(My^, M^^, 0) is very easy 

to calculate using Eg. (13. 53 p . 
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Using fp = 0.152 GeV^ given in [TO], we get: 



0.3 



^ T(15) TT + TT-) 

Similarly 

r(n ^ r,.7) 



^ T(15) TT+TT- 

For Y5 — i- X57 and — )• rj^'j, the lack of experimental results do not allow us to draw any 
solid conclusions. However, we have made their theoretical decay rate predictions. 
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Chapter 5 



Conclusion 



We have analyzed the spectroscopy of particles with hidden charm and bottom in diquark- 
antidiquark structures of the kind of QQ = [qq][qq], where q = u or d ; q = c or b, and Q is 
a diquark. The idea that the color diquark is handled as a constituent building block is at 
the core of the approach taken in this dissertation. In this context we used the constituent 
quark model (CQM) as a benchmark to identify and compare the properties of the newly 
discovered. Constituent quark models typically assume a QCD motivated potential that 
includes a Coulomb-like one-gluon-exchange potential at small separation and a linearly 
confining potential at large separation. In the constituent quark model hadron masses are 
described by an effective Hamiltonian: = + ^2/Cjj(Sj ■ Sj) that takes as input 

i i<j 

the constituent quark masses and the spin-spin couplings between quarks. The coefficient 
K-ij depends on the flavor of the constituents i, j and on the particular color state of the 
pair. By extending this approach to diquark-antidiquark bound states it is possible to 
predict tetraquark mass spectra. The mass spectrum of tetraquarks [ggjfgg] with q = u, d 
and q = c, b neutral states are described in terms of the constituent diquark masses, mg, 
spin-spin interactions inside the single diquark, spin-spin interaction between quark and 
antiquark belonging to two diquarks, spin-orbit, and purely orbital term [6]. We calculated 
the couplings for color singlet combinations from the known Lqq = mesons. Spin-spin 
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coupling for quark-quark in color 3 (antitriplet) state are calculated from the known L = 
baryons. According to one gluon exchange, from Eq. fl2.15p and Eq.f l2.16|) . we have {ICij)^ = 
|(/Cij)o. This relation holds for singlet and antitriplet states. The couplings corresponding 
to the spin-spin interactions have been calculated for the color singlet and color antitriplet 
only. The couplings are not necessarily in the singlet state but octet couplings {}Ccc)g are 
also possible. The quantities JCgg, Kcq and /Ccc involve both color singlet and color octet 
couplings between the quarks and antiquraks in a QQ system (A quark in the diquark Q 
could have a color octet spin-spin interaction with an antiquark in the antidiquark Q). We 
have calculated a relationship: /Ccc ([cq'][cg]) = \ {ICcc)q, using one gluon exchange model. The 
states are classified in terms of the diquark and antidiquark spin, Sq and Sq, total angular 
momentum J, parity P and charge conjugation, C. We have considered both good {Sq = 0) 
and bad {Sq = 1) diquraks for which Lqq = 0, we have six possible states. By diagonalizing 
the Hamiltonian given in Eq.f l3.39|) and using the spin couplings we have calculated the 
masses of the tetraquark states. The mass of the diquark [cq] was fixed by using the mass 
of X(3872) as input, yielding m^cq] = 1.933 GeV. The 1+^ state is a good candidate to 
explain the properties of X(3872). In order to reduce the experimental information needed 
we estimate the remaining diquark masses by substituting the costituent quark forming the 
diquark. Isospin-breaking introduces a mass splitting and the mass eigenstates called ^[c,/] 
and X[c,h] (for lighter and heavier of the two) become linear combinations of X[cu] and X[cd] ■ 
One can put: X[c,i] = cos 9 X^^cu] + sin6' X[cd] and X[c,h] = — sin6' X[cu] + cos6' X^cdi- Besides 
these two neutral states, two charged states arise as a natural prediction of the tetraquark 
picture X+ = [cm] [erf] and X~ = [cd][cu]. The charged partners are X^ are not observed 

We modified the formulism of diquark-antiquark model to calculate the spectrum of 
hidden bottom (bottomness = 0) states for Lqq = 1. We have shown that Yf, is J^'" = 1 
state, with Y[6g] = {[bq]s=o\bq]s=o)p-wa.ve, with the value My^^^ (for q = u,d) equal to 10890 
MeV. We identify this with the mass of the Yb from Belle [13], apart from the T{5S) and 
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T{6S) resonances. 

We discussed the decays modes of X(3872) and y( 10890) states on the basis of quark 
rearangement in the QQ system. Originally the X(3872) was found through its decay into 
J/ipTT^n^ . The decay of a diquark-antidiquark bound state into a pair of mesons can occur 
through the exchange of a quark and an antiquark belonging respectively to the diquark 
and the antidiquark. The only available ones are J/ip — ^ 27r and J/ip ^ Stt, dominated 
by p° and experimentaly confirmed respectively. We have calculated hadronic decay 
widths: r(X[c,i] ^ V + tt+tt-) = 3.78 MeV and r(X[,j] ^ V^ + 7r+7r-7r°) = 0.66 MeV. We got 
information on the mixing angle from the decay rates: ( ) = ('^"'^^ '^^^^^1^ ■ = 0.802. 

° ° \^r(27r) J (cos & —sinOy {p}p 

[c,l] 

But ^f||fjj = 0.8±0.3stat^0.1syst- Thus 9 = ±20° , for X[c,/] and X[c,h] respectively. For 

the charged state X^, that decay via p-exchange only: r(X^ — )■ J/i/jtc^tt'^) = 2\A'^\ J^l^i = 

4.6 MeV. The amplitude for the radiative decay X — )■ J/ip^y is calculated using the Vector 

Meson Dominance (VMD). The radiative decay proceeds through the hadronic transition 

X — 7- Jf-ipp- We found that the ratio of radiative decay to hadronic ones i-e, r(j{^'^/'!/TT+l-) ~ 

0-44 and p^ff^J^g^-) ~ 0-89 ^^at are in agreement with experimental values [28]. We 

exploited the result obtained for the width of X — J/tp'j to give an estimate of the decay 

width into J/ijj'yy. We have obtained r{x^j/4nT+-K-) ~ 3 x 10^^ that is greater than the upper 

limit provided in experimental data |19]. The inconsistency of the theoretical prediction with 

respect to data is not dramatic if one we take into account the very strong assumptions made 

to derive the decay widths. 

For bottomonium systems, the corresponding decay widths are determined by the wave 

functions at the origin for the T{nS), \E'{,5(0), and by the derivative of these functions at the 

origin, \l/'j^g(0), for the P-waves. Due to the possibly larger hadronic size of the tetraquarks 

compared to that of the bb mesons, we modified the Quarkonia potential. For example, 

the Buchmiiller-Tye QQ potential [H]. This will modify the tetraquark wave functions 

^qq(O) from the corresponding wave functions of the bound bb systems, effecting the decay 

amplitudes and hence all the decay widths of the tetraquarks. The corresponding value for 
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the tetraquark states is then calculated taking into account the ratio of the string 

tensions k = ^qq/^^qq, string tension in a diquark fcgg is expected to be different than the 
corresponding string tension kqq in the QQ mesons. We expect k, to have a value in the range 
K E [|;"^]- As the linear part of the confining potential determines essentially the heavy 
Quarkonia wave functions, we find that to a good approximation: \1'qq(0) ~ K\E'ft5(0). The 
partial electronic widths for the P-states are given by the well known Van Royen-Weisskopf 
formula. For the yb(10890) state we have obtained Fee ~ 0.12KeV. This value is close to the 
experimental value given in [31]. We have also calculated the two-body hadronic decays of 
the n(10890), Yb{q) Bj*\k) Bj*\l) . The amplitude for the radiative decay n ^ Xbl will 
proceed through the annihilation of a pair of light quarks into a photon. We have calculated 

r{Yb^-Sr(is)-l+n~) ~ ^'"^ r{Yt^T'(is)t^+TT- ) ~ '^^ which we do not have any experimental 
results. We expect that our predictions should provide guidance for the future. 

This work will be of interest both for theoretical and experimental particle physicists 
for the next couple of years. Now the next task is to explore the tetraquark state [&q'][cg]. 
The motivation is in fact that the searches for tetraquarks in the Be system may find these 
[6g][cg] tetraquark below the BD threshold. Possible exotic signature include strong or 
electromegnetic decays in to B^n or _Bc7, weak decays producing additional peaks in the 
mass spectrum of B^ decay final state. 
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